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A quantized form of Surface Plasmon Polariton (SPP) modes propagating on the metal thin film 
is provided, which is based on the Green’s tensor method. Since the media will be considered lossy 
and dispersive, the amplification and attenuation of the SPP modes in various dielectric media, by 
applying different field frequencies, can be studied. We will also illustrate the difference between 
behavior of coherent and squeezed SPP modes in the amplifying media. 
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I. INTRODUCTION 

The study on surface plasmon polariton [T] is an growing area which has attracted much interest for various 
applications [5]. Because of the quantum nature of SPP EHg, its applications in some areas such as quantum 
information process becomes an active field [7H9]. In order to applying the quantum plasmonic, a suitable quantized 
form of SPP must be provided. In Ref. nni a quantum mechanical form of SPP’s field vectors based on Hopfield theory 
presented but in this formalism the dissipation is not considered. Recently in m for SPP propagating in the semi 
infinite geometry, we have proposed another method for quantization based on Green’s tensor method [EHni which 
contains the loss. Moreover, this method has the potential for generalization to the dispersive and inhomogeneous 
media with different geometries. 

In the present contribution we extend the technique developed in m for quantization of SPP mode for a thin film. 
It also enable us to studying the SPP phenomena in quantum approach such as amplification or attenuation of the 
SPP modes for some quantum (coherent and squeezed SPP) states. 

This paper is structured as follows: the main foundations of quantization of EM fields is provided in section 2. 
By considering the slab geometry, the procedure of evaluating of the corresponding Green tensor and applying it for 
obtaining the quantization form of SPP field vectors are presented in section 3. In this section we also investigate some 
well known relations such as field fluctuations, canonical commutation relations and Langevin equation. In section 4, 
by applying the quantized form of SPP field, we investigate the influence of the frequency and dielectric media with 
different optical parameters on the amplification and attenuation of SPP modes. Moreover, we illustrate schematically 
the difference between behavior of the two modes of SPP( symmetric and antisymmetric) under considered conditions. 
We also compare the behavior of two modes for two types of states that is possible only in the quantum scheme. A 
conclusion is given in section 5. 


II. PRELIMINARIES 

The fundamental concepts that form the basis of the quantization procedure by Greens tensor method are discussed 
comprehensively in [Mm, so a short review is presented in this section. 

The EM-wave propagates in a dielectric medium with dielectric function e(r, w), which is related to the complex 
refractive index 


e(r, uj) = [n{r, w)]^ = uj) + iK{r, w))^, (1) 

here ri{r,uj) and K(r,ui) are real and imaginary part of refractive index, respectively. In general, in a range of 
frequencies, if K(r, uj ) is negative, the dielectric media attenuates the EM- waves, otherwise it can be considered as 
amplifying media. In order to investigating the behavior of EM- waves propagating in the dielectric media, in the 
quantum approach, it is useful to consider the electric and magnetic operators {E{r,uj), B[r,uj)) according to the 
vector potential operators A{r,uj) 


E{r,uj) = 


dA{r,uj) 
at ^ 


B{r,uj) = V X A{r,uj). (2) 

On the other hand, in the frequency domain, the field operators can be considered as positive and negative components: 

E{r,uj)=E+{r,uj)+E-{r,uj), (3) 


where 


- 1 

E^{r,t) = / dujE^{r,uj) exp{jfiu}t). (4) 

v27r Jo 

Accordingly for B(r,uj) and A(r, w) we have the similar relations. By substituting the Eqs. (§ and <§ into the 
quantized Maxwell equation [Biin], and decompose it into the positive and negative parts, a general equation for 
vector potential operator will be obtained, 


- V X V X A+{r,uj) + ^e(r,a;)A+(r,a;) = -noj^{r,uj), 


( 5 ) 
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where j'^{r,oj) is the noise current operator associated with the noise sources in the absorbing (or dissipative) media, 
which is deduced from the fluctuation-dissipation theorem [ISIEI]. It is convenient to express according to 

the normalized noise operator f{r,ui) 

= VQ)(w)/(r,w), 

where the coefficient a{uj) depends on the optical properties of the media, which satisfy the commutation relations, 

,uj')] = 6{r - r')S{uj - uj'), 

[f(.r,u})J{r',u}')] = [P{r,u}),P{r' = 0 . ( 6 ) 

One of the solution of Eq.([^ is based on the standard Green’s tensor method 

/ -l-oo ^ / » / 

dr'G{r,r',a;)J)^(r',w). (7) 

-OO 

The Green’s tensor must satisfy the Eq. (§ when the source is replaced by a point source. 


— V X V X G{r,r ,uj) H— ^e{r,uj)G{r,r ,uj) = IS{r — r ), 


( 8 ) 


where I is a unit tensor. A suitable way for obtaining the system’s Green’s tensor is eigenmodes expansion method 
[HHin] . Considering the eigenmodes and eigenvalues form of Eq. ([^ 


w 


- V X V X Anir,ui) + — e(r,w)A„(r,a;) = e(r, a;)A„A„(r, w). 


(9) 


where A„ and A„(r, w) are eigenvalues and eigenmodes, respectively. The eigenmodes satisfy the orthognality condi¬ 
tion: 


/ -l-oo 

e(r,uj)An(r,uj) ■ [A^(r, a;)]*d^r = 

-OO 


Therefore the Green’s tensor is given by. 


G(r,r') = 


An{r)[An{r’)]* 

NnA„ 


( 10 ) 


( 11 ) 


III. QUANTIZATION FOR SPP FIELD IN A METAL SLAB 

In this section we apply the quantization procedure for a metal slab embedded between two dielectrics, as is shown 
schematically in Figj^ For this system, the dielectric constant can be considered as. 


t ^ 


1 

dielectric 

II 

X 

metal 

1 dielectric 


FIG. 1. Schematic representation of the metal slab embedded between two dielectrics, SPP propagates on each interfaces of 
I-II and II-I. 


e{r,uj) =ed{uj)Q{-z) + 

em(w)0(2:)0(-(z - d)) + ed{uj)0{z - d). (12) 

In order to obtain the vector potential operator (Eq.([^) and quantize it, the Green’s tensor would be derived at the 
first step. 
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A. construction of Green’s tensor for a slab 


By applying the mode expansion method explained in previous section, the Green’s tensor can be obtained. By 
solving the Eq. 0 for the geometry depicted in Fi^ one can find two vector potential modes which correspond to 
two types of SPP modes. They are symmetric and antisymetric modes with the even (lower sign) and odd(upper 
sign) vector potential function, respectively, 


Ak^{x,z) = (x-i— 
^'o 

Ak^{x,z) = 

h 


A{{x + z —z —x e 

V, 


z <0 
0 < z < d 

ikj. 


{x, z) = + z — z>d 

Vo 


where 


(13) 


A = (1 =F exp(-z/md))’ 

2 _ 2 _ 

Vo — ^2 ’ 


= ki- 


(14) 


The SPP modes propagate along the x axis and z/q and Um are the decay coeficients along the z axis for dielectric and 
metal region, respectively. The upper sign in Eq.( 13) is related to antisymmetric mode and the lower sign corresponds 


to the symmetric mode. The frequency of the antisymmetric mode is higher than the frequency of the SPP for a single 
interface while the symmetric mode’s frequency is lower [T] . Furthermore, odd and even mode’s frequency satisfy the 
dispersion relation: 


^ -r , 


- 1 


^dk'm 


+ 1 


(15) 


On the other hand, by inserting the Eq.( 131 in Eqs.Q and (101, the eigenvalues A„ and the normalization coefficients 
Nn can be obtained. 


7V„(fc,)=27r{^(l + %) + 
Vq v. 


0 

A ^^{1 + %)(1 - ^ 


Vm Vt, 

=2t: N^{k^) 


^)T2d(l-^)e- 


\ —Ip' 

An — 


(16) 


According to Eqs. (13), (16) and (11), the Green’s tensor can be written as: 

^ikx{x—x ) 


G{r,r',uj) = x 


j^2 _ , 2 

^ + z0+ 


X 


{{x - i —z)e"»"0(-z) T (i + Z— - d) + 
Vo Vq 


A[{x + i — z)e 

Vm. 


=F (x — z — z]e 


kx 

Vr, 


_,4^$)e-™G-d)]0(z)0(d-z)}x 


{(x - i — z)e''°^'Q{-z) T [x + z —z)e-"“("'- 
Vo Vq 


- d) + 


kx 

Vrr 


A[{x + i — z)e =F 


(x - z —z)e"’"("'-'^)]0(z')0(d - z')} 


(17) 
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In order to evaluate the integral by residu theorem, the poles of the denominator must be evaluated. For analytic 
solution, it is convenient to consider the kx very close to the roots of the denominator. We assume the are the 
roots, the upper (lower) sign is correspond to odd (even) mode. By taylor expansion for Vm, about the roots, the first 
order approximation is given by m- 


l^mikx) — Vmik^pp) + {kx k^pp) \k^ = kfp, 


By applying Eq.(18), the denaminator of Eq. 0 yields: 


Em^O “f ^m{kx) k^pp — i^x k^pp) ^ { i^x + ^Jp)d“ 


dkx 


(18) 


(19) 


As can be seen, the poles of the denominator are the zeros indeed. Eor simplicity in the subsequent calculations, we 
use kf instead of kf So by some algebra, the Green’s tensor can be obtained as: 


G(r,r',w) = 


{(i - i’^z)e''°^e{-z) T (i + i'^z)e-'°^^-‘^'>e{z - d)+ 

Vo Vo 


.kf 


.kf 


A[{x + i'-^z)e-'’-^ T (5 - z —z)e^™(^-‘^)]0(z)0(d - z)}x 

Ur) 


{{x - e{-z ) T (i + i — -‘^^e{z - d) + 

Vo Vo 

A[{x + i^z)e-'^^' t{x- z^z)e"’"("'-‘^)]0(z')0(d - z')} 

Urn. Urn 


( 20 ) 


where 


D = 


2{(-K(fcf)(i^m(A:f)^|fcT +fcf)} 


( 21 ) 


Also the z^o’s and z^m’s in Eq.(20) are evaluated for kfpp. 
potential operator on a slab. 


Now by using the Green’s tensor, we can quntize the vector 


B. Field quantization and canonical commntation relation 


Since, in the system at hand, there is three regions (dielectric- metal- dielectric), the corresponding 3 components 
of noise current are given by, 

j^{x,z,u:) =j^+(a;,z,a;)0(-z) -f j^+(a;,z,a;)0(z - d) + 
j]^+(x, z, uj)Q{z)Q{d - z), 

= [\l a‘^{uj)Q{—z) + a’^{ijj)Q{z — d)+ 
y/a"^{uj)Q{z)Qid - z)]f{x, z,uj). ( 22 ) 


By applying the Eqs. (22), (20) and ^ and some calculations, one can obtain the vector potential operator and 
introduce the annihilation and creation operators like the procedure in ?? 


i+(r,a;) = inoD{^^^)^ x 

{{x - z^z)e"“"0(-z) T (i + i —^)e-"“("-'^)0(z - d)+ 

Vo Vo 

A[{x + z^z)e-"-" T (i - i —z)e"™("-'^)]0(z)0(d - z)}x 
{a^{x,uj) +dj^ix, uj)}, 


(23) 
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here is the imaginary (real) part of the and 


Piu;) =|aV)l(l + ^) 


|a™(a;)P|^x{(l 

^ I., 12'' 


|fc?l\ 2 
'^1^0 + 1^0 

|fcT|2 2(l-e-("-+‘'-)^) 


— ut^d _ ^—^md 


}■ 


(24) 


^rn 


In Eq.(23l the operators and dj indicate the annihilation of the SPP mode with symmetric (+) or antisymmetric 
(—) field function that propagates rightwards and leftwards respectively and have the explicit form as: 

olT / 

aT(x,cc) = (^)U / dx'dz'e*'=f("-")j+(a:',z',a;)- 

{(i - i^z)e''°^'Q{-z) T (x + - d) + 

Vo Vo 


A[{x + P^z)e 

^m. 


T 


{x - z^z)e""*("'-'^)] 0 (z') 0 (d - z')} 


(25) 


and 


df{x,Lo) = {^)-^ f f dx'dze-^’^f^--^'^jUx,z,ujy 

P H J-oo Jx 

{(i - z^z)e^“^'0(-z') T (x + i—z)e-'°^^'-'^^eiz - d)+ 

Vo Vo 

A[{x + z^z)e-"-" Tix- z —-‘^)]0(z )0(d - z )} 

^m. ^m. 


(26) 


By some calculations one can find that the annihilation and creation operators satify the commutation relation, 

[a^{x,u}),Ap{x ,uj)] = [dy{x ,uj ),al\x,uj)] = 

5(uj — uj ) exp(zA:^(x — a; )) exp(—fcjja: — x |), (27) 

[a^{x,u}),dmx ,uj )] = [d^ix ,uj ), a^'’'(a;, w)] = 


f ' 2k^ ' ' 

6 {u} — oj )Q{x — X )^p- exp{—kfj{x — x )) sinkfj^{x — x ). 

^sR 

On the other hand, one can explore the canonical commutation relation and obtain, 

[A{r,t),-eoE{r',t)]= f dw.— ^ ■ • ImG’(r,r',a;), 

Jo TTCoC^UJJ (w) 


(28) 


(29) 


The details of the calculations and the explicite form of 7 are given in Appendix A. By considering the general 
property of Green’s tensor lim|,j|_,.oo ^G(r, r , 0 ;) = —S{r — r ) and assuming that: 


/3 (w) = 2/ieow^7 (w) 

one can prove that the canonical commutation relation is satisfied. 

[A(r, t), -eo£^(r , t)] = ihS{r - r ). 


(30) 


(31) 


On the other hand by substituting Eqs.(24) and (41) into Eq.(30) the explicit form of a”^{uj) and a‘^{uj) can be 
obtained. 


| q ;™( w )| = 2 huj^ eolmem{i^), 

\a'^{uj)\ = 2da;^eoImed(a;). (32) 
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IV. NUMERICAL RESULTS: EXPLORING THE AMPLIFYED AND ATTENUATED CONDITIONS 
FOR COHERENT AND SQUEEZED SYMMETRIC AND ANTISYMMETRIC SPP MODES 

A. magnetic field 


In order to investigate the variation of the SPP’s modes propagating in amplifying and attenuating media, we need 
to calculate the magnetic field. For simplicity, we consider that the SPP waves propagate rightwards. By applying 
Eqs. (25), (23) and (§ magnetic field operator can be obtained: 


H+{r,uj) = a^(x,w)x 

LT2 LT2 

{(i^o - ^)e"«^0(-z) T (-1^0 + - d)+ 

LT2 

A[i-,ym + T Wm - - z)}, 

^m. ^m. 


(33) 


where a^(a:,w) satisfy quantum Langevien equation: 


dd^{x, oj) 
dx 


= ikf a^{x,uj) + F'^{x,uj), 


(34) 


where 


Oj^T /-oo 

F^(x,w) = (^77^)5 / dz j(^(a;,z ,0;)- 

P [^) J-oo 

{(i - f^z)e"«"'0(-z') Tix + i^z)e-''A^'-‘^'>Q{z - d) + 
t'o ^'o 

A[{x + t{x- i^z)e"’"("'-'^)]0(z')0(d - z')} (35) 


here F^{x,uj) is the operator associated to the Langevien noise source. On the other hand by attention to the explicit 
form of F^{x,uj) one can find that this operator can only appear in the absorbing(or dissipative) media when noise 
sources are exist. 

By solving the Eq.(34| and considering the general property of noise source operators (j(^(a;, z, w)) = 0, the average 
form of Eq. (33) is given by 


(iL+(r,a.)) = x e*'=?-(a^(a;))x 


LT2 LT2 

{(i^o - ^)e"«^0(-z) T i-i^o + ^)e-"“(^-‘^)0(z - d)+ 

r'o 

j!,=F2 ^T2 

A[i-,ym + T Wm - ^)e''-^^-‘^^]Q{z)Q{d - z)}. 


(36) 


The average of magnetic field operator can be considered for the different kinds of SPP states like coherent and 
squeezed states [51 E]. The annihilation and creation operators of the SPP modes obey the relations of bosonic 
operators. Therefore, when the SPP states are prepared in coherent and squeezed states the following relations can 
be considered. 


afl(w)|a) = a\a). 

aR{oj)\^, a) = - i/a(a;))|a). 


(37) 

(38) 


where a = |a|e*®, fi = cosh(|^|) and ly = sinh(|^|)e^^^. |^| and 6^ are the absolute value and argument of squeezed 
parameter, respectively. By applying the Eqs.(|37| and (38), the magnetic held average can be calculated for coherent 
and squeezed SPP states. 














B. Studying the influence of frequency on the SPP modes 


The optical properties of the dielectric media adjusted to the metal film can affect the properties of SPPs. Ac¬ 
cording to the intrinsic absorption property of the metal, SPP modes suffer damping which reduces the SPP length 
propagation. On the other hand the dielectric media with negative imaginary part of refractive index {n = Ud + ikd) 
can act as gain media. When the gain of dielectric media is sufficient to compensate the loss in the metals, the 
will be negative and the system acts as amplifying media otherwise for positive k^j the system can be considered as 
attenuated media [2H1 - EH] . 

The variation of k^j according to the frequency can be shown by the dispersion relation diagram. It is depicted in 
figureIn Fig. [^despite of different behavior for symmetric and antisymmetric modes, there is a same prediction. 



FIG. 2. Symmetric and antisymmetric SPP modes dispersion relation, (a) real part and (b) imaginary part of the wave 

number. Where d=60 nm, the dielectric media with n = 1.9726 — iO.081 and the metal film with ujp = 14.02 x 10^®(Rad/s) 

and 7 = 6.25 x 10^®(Rad/s) have been considered. (The data are in [22]). 

In general, in the frequency ranges that in the loss case kfj is negative, the SPP modes can be amplified otherwise 

they are attenuated. 

Moreover, it can be shown that for thick film the dispersion relation of symmetric and antisymmetric modes are 
identical and are in accordance with the one interface SPP mode [23] . 


C. Studying the influence of the dielectric media on the SPP modes 

Besides the frequency, the different gain media (dielectric media) can affect the SPP modes behavior. This is 
illustrated schematically in Fig. 

Fig|^ shows some interesting points: hrst, any gain media adjusted to the metal thin film can not amplify the SPP 
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FIG. 3. Variation of ksi of the SPP modes in the vicinity of different gain media, where Ud = 1.9726 and uj = 4.8 x 10^^Rad/s. 

modes. Indeed, for amplifying kd and kfj must be negative simultaneously. Second, for a particular kd the operation 
of the system is different for symmetric and antisymmetric modes. This is discussed in more detail later. Third, the 
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slope of the lines indicates that for a given kd range, the variation of the symmetric mode in comparison with = 0 
(no gain and no loss) is very greater than the antisymmetric mode. Fourth, except of very small range of kd (very 
close to the = 0 for symmetric mode) the rate of amplifying or attenuating symmetric mode is very faster than 
the antisymmetric mode. 

In order of illustration and comparison the variation of SPP modes in these ranges, we plot the magnetic field average 
Eq. (36) for coherent symmetric and antisymmetric modes for different ranges of kd ■ 


1. Investigating the variation of coherent SPP modes under different gain media 

According to Fig[^ for different ranges of kd the SPP modes suffer different conditions. The first condition is 
where two SPP modes are attenuated. For instance, by considering the dielectric media with n = 0.9726 — i0.063 
the k'fj of the SPP modes are positive. It means that the gain of the dielectric media can not compensate the loss 
of the metal and the SPP modes are also attenuated. It is shown schematically in Fig|^ . In the second condition. 



FIG. 4. Attenuation of two modes on the z = 0 interface, where \a\'^ = 7,6 = 1.5 Rad. fcj} = 7.8029 x 10^(Rad/nm) and 
k- = 2.7617 X 10®(Rad/nm) . 

the system’s operation is different for symmetric and antisymmetric modes. For example, the dielectric media with 
n = 0.9726 — t0.072 cause the system attenuates the symmetric mode but amplifies the antisymmetric mode. It is 
shown in Figj^ The amplification of symmetric and antisymmetric modes simultaneously is occurred in the third 



FIG. 5. Attenuation of symmetric modefcj} = 2.5858 x 10^(Rad/nm) and amplification of antisymmetric modek^j = —1.96 x 10® 
(Rad/nm) on the z = 0 interface. 


condition where for a given kd, the is negative for two modes. Fig. shows this condition where the refractive 
index of dielectric media is n = 0.9726 — 10.08. It means that for the two modes the gain of the dielectric media is 
sufficient to overcome the loss of the metal. 

As the Fig [^illustrates except the small range of kd,the magnitude of fcj} is greater than kjj, therefore for a special 
condition the rate of variation for symmetric mode is noticeable than the antisymetric mode, as Figs.([4]|^ demonstrate 
obviously. Moreover by comparison the Figs.(4][6 1, one can find that for different conditions, the changing of symmetric 
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FIG. 6. Amplification of two modes on the 2 = 0 interface, where = —1.8673 x 10^(Rad/nm) and = —6.1217 x 10® 
(Rad/nm). 


mode behavior is more noticeable than another. 

A very interesting point is that, only for one magnitude of the behavior of two modes is the same. It is the 
point of intersection of two graphs in Fig|^ Figj^ shows this condition where the optical property of dielectric media 
is n = 0.9726 — 10.07747. Moreover, the same figures are obtained for the modes that propagate along the z = d 



FIG. 7. The same behavior of two modes on the 2 = 0 interface, where fcj} = k^j = —4.81 x 10®(Rad/nm) . 


interface. 


2. Comparison between coherent and squeezed SPP modes state 


By applying Eqs.(37) and (38) into the Eq.(36), the average of magnetic field is obtained for coherent and squeezed 
SPP state. For more investigation, we consider a = |a|e®® as a complex number. The influence of the phase 6 on the 
average of the SPP magnetic filed for one interface system is studied in m- Such result is obtained for two interfaces 
system which shows that the difference between coherent and squeezed state is significant for 0 = 1.5. In Figj^ 
the difference between the magnetic field average for coherent and squeezed states that propagate in the amplifying 
system is illustrated. The average of magnetic field for antisymmetric mode is shown in Fig|^ From Figs. and 
one can deduce that the drastic difference is occurred between squeezed and coherent states. It can be generalized to 
other amplifying or attenuating system. 


V. CONCLUSION 

In this paper we have provided another approach for quantization of SPP modes on the thin film structure, based 
on Green’s tensor method which contains some new quantum concepts in the SPP field such as noise current, field 
fluctuations, and Langevin equations. Moreover, this approach enable us to study the influence of the different 
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FIG. 8. The difference between the average of the symmetric mode’s magnetic field, where n = 0.9726 — i0.08 is chosen for 
dielectric media. 



FIG. 9. The difference between the average of the antisymmetric mode’s magnetic field. 


conditions on the propagation of the SPP modes. The results are classified as follows: 

First, for certain media, the variation of frequency can cause SPP modes amplification or attenuation. It also has 
been shown that the behavior of symmetric and antisymmetric modes is different in frequency domain. 

Second, for certain frequency, the amplifying or attenuating SPP modes is dependent on the optical parameter of the 
dielectric media adjusted to the metal film. We have also compared the behavior of two SPP modes with each other 
for different media. 

Third, we have illustrated that the drastic difference is between different types of SPP modes, i. e., coherent and 
squeezed states. 
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APPENDIX A 


In order to prove the Eq.(29), it is neccesary to derive some relations: 


2 (w) 


[A{r,t),-eoE{r ,t)]= duj\D\ ' x 


^sl 


|jL=F|2 ikf\x—x I —ik*^\x—x \ 

}x 




kf 


fvg 


{(i - ^^5)e"'>"0(-z) Tix + - d)+ 

I'D t'O 

A[(i + T (i - i—i)e"'"("-'^)]0(z)0(d - z)}x 

{(i + z^z)e"o"'0(-z') T (i - i^z)e-''o(^'-d)Q(z' - d) + 

^0 ^0 

A*Ux - T (i + i —z)e"-("'-‘^)]0(z')0(d - z')}. 

u* i/* 


(39) 


On the other hand, the Green’s tensor for a dielectric-metal- dielectric structure (see Eq.([20|)) satisfy the following 
relation: 


dslme(s, w)G(r, s, uj) ■ G*{s, r ,uj) = 


here 


\kf\^ 

^kfnkJi 


7 {u;)\Dm- 


„ik^\x—x I p — ik*^\x—x \ 


kj 




-}x 


{(i - i^2)e"'>"0(-z) Tix + - d)+ 

A[{x + Tix- i^5)e'^”(^-‘^)]0(z)0(d - z)}x 

{{x + i^z)e‘'o^'ei-z) t{x- i^z)e-''oG'-d)Q(z' - d) + 
A*\{x - i^z)e-'^-^' T (i + z^z)e'^™(^'-‘’*)]0(z')0(d - z')}, 

V* u 

m xy-i 


7 (w) =Imer„(l-f 


1^0 + ^^5 


ImedlAI-^ X {(1 + 


|yfcT|2 2(l-e-("™+"-)‘^) 


(40) 


T(l-^) 


|fcT|2 2(e-'^™‘'-e-'^-'^) 


}• 


(41) 


ande(s,CLi) has been given in Eq.( 12). By rewriting Eq.([3^ according to Eq.(40) and considering the general property 
of Green’s tensorlMl 


dslme(s, a;)G(r, s, w) • G*(s,r ,oj) = -^IniG(r, r ,uj) 


(42) 


The desired equation can be deduced: 


[i(r,t),-eoE(r ,t)] = f -^ImG(r, r ,w). 

J Try (w) UJ^ 
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